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The electronic structures of graphene systems and topological insulators have closely-related fea- 
tures, such as quantized Berry phase and zero-energy edge states. The reason for these analogies is 
that in both systems there are two relevant orbital bands, which generate the pseudo-spin degree 
of freedom, and, less obviously, there is a correspondence between the valley degree of freedom in 
graphene and electron spin in topological insulators. Despite the similarities, there are also sev- 
eral important distinctions, both for the bulk topological properties and for their implications for 
the edge states - primarily due to the fundamental difference between valley and spin. In view of 
their peculiar band structure features, gapped graphene systems should be properly characterized 
as marginal topological insulators, distinct from either the trivial insulators or the true topological 
insulators. 



I. INTRODUCTION 

Condensed matter physics is witnessing an increasingly 
rapid development marked with a sequence of surprising 
discoveries. Two good recent examples are the realiza- 
tion of graphene^ ~ a plane of carbon atoms forming hon- 
eycomb lattice, and topological insulatori^MI _ materials 
that are insulating in the bulk but conducting at the 
surface owing to topological reasons. In fact, the two 
examples were linked since the beg inn ing: quickly after 
the initial experiments on graphend^, it was suggested 
by Kane and Mele, that graphene may be gapped due to 
the intrinsic spin-orbit interaction and provide a proto- 
type of a novel class of time- reversal-invariant topological 
insulatorPISl. 

This idea has since inspired numerous works that have 
surprisingly, yet significantly, complemented our knowl- 
edge of some seemingly well developped fields of solid 
state physics. The original prediction by Kane and Mele, 
however, remains unfulfilled in graphene because of the 
very weak spin-orbit interaction in realistic graphene 
samples. 

In fact, in graphene, spin-orbit interaction is so weak 
that in most existing experiments without a magnetic 
field, the spin degree of freedom hardly plays a discernible 
role except for contributing an additional degeneracy. 
This certainly prevents graphene from being an ideal rep- 
resentative of (time-reversal-invariant) topological insu- 
lators for which strong spin-orbit interaction is normally 
a crucial ingredient. Despite this, graphene and topo- 
logical insulators share a number of clo sely related fea- 
tures like a quantized Berry phas^^*^ and zero-energy 
edge stateaiSHSl Both in graphene and topological insu- 
lators one can isolate two relevant bands with the oppo- 
site orbital symmetry, which in the band structure are 
connected by the matrix elements linear in momentum, 
leading to the Dirac cones. Further, in both systems, a 
gap in the Dirac dispersion can be opened, by means of 
spin-orbit interaction in topological insulators, and sim- 
ply by creating a sublattice on-site energy imbalance in 
graphene. The analogy between the two systems is com- 



pleted by associating electron spin in topological insula- 
tors with the valley degree of freedom (also time-reversal 
odd) in graphene. At the formal level, as we shall see 
in a moment, the effective theories for valleys and for 
spins bear such great similarities, that one could very 
well be tempted to transplant the conclusions obtained 
in topological insulators into graphene systems, with spin 
substituted by valley. This raises the question: how far 
exactly can such a spin and valley analogy be pushed? 

To examine carefully the analogies and the differences 
between valley-based graphene-like systems and spin- 
based topological insulators, it is helpful to compare spe- 
cific examples which nevertheless possess generic proper- 
ties. For the graphene case we will use electrostatically- 
biased bilayer graphene (BLG^i^, where an electronic 
bulk band gap can be opened in a practical and con- 
trollable manneJl^^i^. For the topological insulator case, 
we will use the Bernevig, Hughes and Zhang (BHZ) 
modeP^, which is a good prototype of two-dimensional 
(2D) topological insulators and has its realization in 
HgTe quantum wellJ^^^^. Our comparison will be fo- 
cused on the characterization of bulk properties and their 
relation to the presence of subgap edge modes. This 
scheme follows intimately the fundamental logic underly- 
ing the idea of topological insulators, namely, the bulk- 
edge correspondence.f^SlSSI ^pj^g outcome of this investi- 
gation has direct experimental implications in terms of 
subgap edge transport. 

We emphasize that the purpose of this paper is by 
no means to provide a stringent and comprehensive re- 
view on the comparison between graphene systems and 
topological insulators - as one of the limitations, we will 
ignore in the following the spin degree of freedom in 
graphene. We highlight several points that illustrate the 
analogies and the differences behind valley numbers and 
spin. The manner of our presentation is in tended to be 
heuristic. We refer the readers to Ref.'^^'^ for more de- 
tails and more rigor. 



II. ANALOGY BETWEEN GRAPHENE 
SYSTEMS AND TOPOLOGICAL INSULATORS 

In this section we demonstrate the analogy between 
graphene systems and 2D topological insulators by com- 
paring two examples: gapped BLG and the BHZ topo- 
logical insulator. In both examples, we consider for the 
moment the ideal cases - meaning, for gapped BLG we 
ignore coupling between valleys, and for the BHZ topo- 
logical insulator we assume spin is a good quantum num- 
ber. 



A. Effective Hamiltonians 

We start by writing down the low-energy effective 
Hamiltonian for each system. For gapped BLG, this 
reads 



HsLcik) = 
Hxik) = - 



Hiiik) 
HK'{k] 

{K+ikyf -A 
HK'{k) = H*K{-k), 



(1) 

(2) 
(3) 



where K and K' denote the two valleys, k is the wave vec- 
tor relative to each valley, and A is the electrostatically 
tunable gap; for the BHZ topological insulator, it reads 
(with immaterial simplifications of the original model) 



HBHz[k) = 
Hs{k) = 



Hsik) 
Hs'ik) 



a{kx + iky) — (m — k^) 
Hs'{k)=H*s{-k), 



(4) 

(5) 
(6) 



where S and S' denote the two spin sectors, to > and 
a 7^ are two real parameters. Note that, for simplicity, 
both Hamiltonians above have been made dimensionless 
by choosing fixed length and energy scales in each model. 

One immediate observation can be made regarding the 
similar roles played by valley and spin. In the BHZ topo- 
logical insulator, time reversal symmetry is respected and 
imposes a constraint on the two spin sectors given by Eq. 

The same constraint, given by Eq. ([S]), governs the 
two valley sectors of gapped BLG also as a consequence 
of time reversal invariance (as mentioned above, due to 
lack of any significant spin-orbit interaction spin degree 
of freedom is decoupled from the orbital ones; therefore 
one can treat electrons as spinless particles in all calcula- 
tions, reinstating the spin at the end). These constraints, 
generally obeyed in the two systems^, have crucial im- 
plications for the overall bulk and edge properties as we 
will see in due course. 

The two effective Hamiltonians, when characterized 
using the usual topological measures, seem to show an 
even closer analogy. For this purpose we use the com- 
mon expression for the topological invariant associated 




FIG. 1. The (pseudo-)spin textures for the eigenstates of: a 
trivial two-band system (top panels), one spin sector of the 
BHZ topological insulator (cental panels), and one valley of 
gapped BLG (bottom panels). The left panels show the band 
structures at low energy and the spin polarization for some 
selected states; the right panels show more details of the lower 
bands. Note that, for illustration purpose, the Hamiltonians 
H{k) — — g(fc)-cr (with an extra minus sign as compared with 
those in the text) are used such that the spin polarization of 
the lower band states is aligned with g. The trivial system is 
defined by setting g(fc) — (0, 0, m + fc^) with m > 0. 



with a fully-gapped two-component Hamiltonian H(k) — 
g(fc) • cr. Here cr is the vector of the Pauli matrices and 
g(fc) is a real vector. In terms of the normalized vector 
g — s/|g| the expression for the topological invariant iJ^S 



(7) 



The number c characterizes the mapping from the 2D k- 
space to the 2D parametric sphere subtended by g (see 
Fig. [ij for details see Section HI I . It is straightforward to 
find for individual valley/spin sectors (assuming A < 0; 
c-KjK' changes sign when A changes sign) 



CK 

cs 



-CK' 

-cs' - 



1. 



(8) 
(9) 



Namely, the signs of c's reverse with respect to val- 
leys/spins, echoing the underlying symmetry; the mag- 
nitudes of c's appear to be the same integer 1 in both 
gapped BLG and the BHZ topological insulator! 
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m > 0, a = 



0<a«m 



m«a 



FIG. 2. Evolution of the band structure corresponding to the 
Hamiltonian ([5| by tuning a, with fixed positive m. This il- 
lustration is to be compared with Fig. [3] to provide a naive 
analogy between the gap opening in BLG and that in the 
BHZ insulator (see text for details). Note that the compari- 
son between a and m should be understood as being under a 
certain fixed length scale which converts the two parameters 
to the same dimension. 



B. Band inversion 

The nontrivial bulk properties characterized by Eqs. 
([s]) and (|9| are closely related to the occurrence of band 
inversion. Here we present a simple account of the band 
inversion in each system, bearing in mind that such ac- 
counts are by no means descriptions of the microscopic 
mechanisms, but rather phenomenological illustrations 
based on the effective models. 

We start with the BHZ topological insulator. Due to 
the symmetry, it is sufficient to look at one spin sector, 
Hs say. Let us first compare two cases, with m < 
and TO > respectively, and a — 0. It is obvious that 
the eigenstates in both cases are completely polarized in 
terms of the two-component pseudo-spin. Meanwhile the 
crucial difference is that in the latter case (to > 0), the 
two parabolic bands intersect each other, while in the for- 
mer case (m < 0) they are fully separated. In a purely 
phenomenological sense, we call such a reversal of the 
order of bands at low energy "band inversion" . But the 
real excitement about band inversion comes only after 
we turn on a, which leads to an avoided crossing be- 
tween the two bands resulting in the opening of a gap 
(which makes the bulk of the system insulating). Fig. [2] 
shows how the profile of the inverted bands changes with 
increasing magnitude of a. Clearly enough, at small k 
(where only to plays an important role), the electronic 
eigenstates tend to preserve their inverted character - 
states in the upper/lower band are largely polarized (in 
terms of pseudo-spin, see also Fig. [T]) the same way as 
the original downward/upward parabolic band; at large 
k (where both to and a become irrelevant), the order 
restores and the eigenstates in the upper/lower band de- 
viate little from the original upward/downward parabolic 
bands; in-between, as long as a ^ 0, the eigenstates con- 
tinuously evolve such that the pseudo-spin polarization is 
gradually reversed. Hence the band inversion has a direct 
consequence on the nontrivial evolution of the pseudo- 




0<t^«\V\ 
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FIG. 3. Evolution of the band structure corresponding to the 
Hamiltonian (10 1 by tuning t±, with fixed nonzero V. This 



evolution illustrates the unusual aspect of the gap opening in 
BLG (see text for details). 



spin polarization that is characterized by the number c. 
More specifically, the occurrence of band inversion leads 
to the occurrence of a topological phase transition with 
cs/S' changing from to ±1. 

In gapped BLG, a similar occurrence of band inver- 
sion can be best seen at the level of the single-valley 
four-component Hamiltonian H^^^^ -corresponding to 
two layers each with two sublattices- from which the 
two-component one in Eq. ^ is derivecP^. Again, the 
symmetry makes it sufficient to consider one of the two 
valleys. Explicitly, we have 



H 



(4x4) 
K 



(fc) = 



V 

hvpk^ 



hvpk- 
V 
t± 



—V hvpk- 
hvFk+ —V 



(10) 



where ±1^ are the electrostatic potential on two graphene 
layers, t±^ is the interlayer nearest-neighbor hopping en- 
ergy (assuming Bernal stacking of BLG), k± — k^ iky, 
vp — 3tao/2h with oq and t being the intralayer nearest- 
neighbor distance and hopping energy. Note that, unlike 
Eq. ([2]), this Hamiltonian is not dimensionless, and the 
parameters here are related to those in Eq. ([2| by defin- 
ing A = V/t± and redefining iok — k with io = hvp/t±. 

We imagine first that the coupling between the two 
graphene layers is turned off {t± = 0) such that, due 
to different electrostatic potentials, the electronic band 
structure at one single valley consists of two shifted Dirac 
cones as illustrated in the left panel of Fig. [Sj For clarity 
we mark with two different colors the initial conduction 
bands (blue) and valence bands (red) in Fig. |3] Then 
we gradually turn on t±. Analogously to the turning 
on of a in the BHZ model, a finite t± opens a gap in 
the bulk and turns the system into an insulator (Fig. 
[3]). As a result, at low k, a small part of the valence 
band of one layer is now "glued" with the conduction 
band of the other layer to compose the new low-energy 
conduction band. The opposite happens to the new low- 
energy valence band (see the middle panel of Fig. |3|. 
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This is seemingly analogous to the band inversion we 
have seen in the previous example. Upon increasing the 
magnitude of the interlayer coupling, the two high-energy 
bands are lifted further apart (the splitting is roughly 
2t± when t± is large); the two low-energy bands each 
approach a parabolic shape, in agreement with Eq. [2] 
while the inverted-band character nevertheless remains 
(see the right panel of Fig. [3| . 

To summarize, the band inversion in gapped BLG, as 
illustrated above, happens as long as ^ 0, and results 
in a finite ck/k' that changes sign ±1 when V changes 
sign. This is similar to the topological phase transition 
exhibited by the BHZ model, but not identical, since in 
the BHZ model cg^gi changes from to 1 as parameters 
are varied, and not from —1 to +1. 



III. LIMITATIONS OF THE ANALOGY 

One fundamental limitation of the analogy between 
graphene and topological insulators, with valley playing 
the role of spin, is that valleys, which correspond to only 
parts of the Brillouin zone, are only defined in the energy 
window much smaller than the electron bandwidth (frac- 
tion of eV). That is in contrast to electron spin which 
is well-defined throughout the Brillouin zone. In most 
cases, due to the symmetry, the two valleys contribute 
oppositely to the topological property of the entire Bril- 
louin zone - Eq. ([8| for gapped BLG is an example. This 
results in overall topologically trivial graphene systems 
(according to the definition of a topological insulator; in 
the present context, this means that the integral in Eq. 
(7) needs to be calculated over the entire Brillouin zone 
and not only around one valley). This is certainly im- 
portant since gapless edge states, if they exist at all, are 
not protected against intervalley scattering. More impor- 
tantly, even in the complete absence of intervalley scat- 
tering (in the bulk and at the edges) crucial differences 
exist between the two systems at the level of single- valley 
models. These differences are the subject of this section. 



A. Characterization of bulk properties 

In the previous section we have employed Eq. ([T]) for 
the two models that we are considering, suggesting that 
the numbers c characterize the topological properties as- 
sociated with bulk bands. We now go back to examine 
the validity of such characterizations, based on a geomet- 
ric interpretation of Eq. ([?]). 

To begin with, we parametrize the general two-band 
Hamiltonian by setting g = (sin 9 cos (p, sin 9 sin (p, cos 9) , 
where 9 and ip are the angles defining the direction of k. 
It is straightforward to show that the integrand of Eq. 
^ is equal to sm9{dk^9dk^(p - dk^(pdk^9). Recognizing 
that the factor in the parenthesis is but the Jacobian 



\d{9,(p)/d{kj:,ky)\, Eq. ^ can be rewritten as 
1 



c = — / / sm 

4tt 



(11) 



It is immediately clear that geometrically c measures the 
solid angle (divided by An) covered by g when k runs 
over the domain of interest. 

Keeping this in mind, we can easily understand, by 
referring back to Fig. [l] that c = ±1 for the BHZ topo- 
logical insulator because g covers the full sphere once, 
while c = ±1 for the gapped BLG because g covers a 
hemisphere (the upper or lower hemisphere, depening on 
the sign of A) exactly twice. This fact reflects the cru- 
cial difference of how the vectors g behave at fc — >■ oo in 
the two models. In the BHZ topological insulator case, 
g(fc — >■ oo) converges identically to the south pole regard- 
less of the direction along which infinity is approached, 
which allows for the compactification of the infinite k- 
plane as a Riemann sphere. In this topologically 
nontrivial mapping from the compactified A:-plane to the 
sphere subtended by g is well defined and c is the topo- 
logical invariant associated with this mapping. In the 
gapped BLG case, however, we have a different behavior: 
g(fc — >■ oo) sits on the equator and varies with the polar 
angle of fc, preventing a proper compactification of the 
fc-plane. Consequently the numbers c obtained , albeit 
nonzero, cannot be identified as topological invariant^^ 
as in the previous case. 

Despite this important difference, an alternative in- 
terpretation of Eq. ([t]) in terms of Berry phase and 
Hall conductivity remains valid, as we proceed to dis- 
cuss. This indicates that the number c has a well-defined 
physical (albeit not topological) meaning associated with 
the properties of the bulk. The Berry phase acquired by 
an electron transported adiabatically along a closed tra- 
jectory C in the momentum space is given by 



J = I ij) {u{k)\V ku{k)) ■ dk 



(12) 



where u{k) is the periodic part of the Bloch wave function 
and can be obtained by taking the eigenstate of the Bloch 
Hamiltonian H{k). Using Stokes' theorem it is easy to 
shown that the above formula is equivalent to 
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7 



5c 



d^kAlm{{dkMf')\9kMk))), (13) 



with Sc the area enclosed by C. For H{k) = g(fc) • cr 
with g(fe) = {sm9 cos (fi, sin 6* sin cos 9) (the same 
parametrization as used above; we consider a fully 
gapped system where |g(fe)| ^ for all fc), the lower band 
eigenstate can be written as u{k) ~ (sin f e^'*^, — cos |)"^. 
Then it is straightforward to show that the integrand in 
Eq. (13) is again equal to sm9\d{9,Lp)/d{k.j;,ky)\, hence 



7 — 2'Kc . 



(14) 



It follows from Eq. ( 13 ) that 7 can also be identified as 
the Hall conductivity(in units of e^/27r/i) contributed 



by all the occupied Bloch state^^. Therefore, Ck/k' 
in Eq. Q, despite not representing a true topologi- 
cal invariant, still corresponds to the valley-specific Hall 
conductivitjl^itt^s]^ (^y^ units of e^//i) when the Fermi en- 
ergy lies in the bulk band gap (i.e., c is the contribu- 
tion of the occupied states in a given valley to the bulk 
Hall conductivity). Similarly, of course, Cg/gi in the BHZ 
model represent the spin-specific Hall conductivity, next 
to being well-defined topological invariants. 



B. Bulk-edge correspondence 

For the BHZ topological insulators cs/s' arc topolog- 
ical invariants. According to the principle of bulk-edge 
correspondence, this implies the presence of one pair of 
spin-helical gapless edge modes at each edg d^^ l '^^ l - their 
helicity follows from the opposite signs of cs and cs' ■ In 
contrast, a similar conclusion cannot be drawn for the 
valley-specific Cx/k'i since these quantities are not topo- 
logical invariants. In general, therefore, we cannot expect 
bulk-edge correspondence in this case. In the particular 
example of gapped BLG, this leads, as we will show, to 
the possibility that the number of gapless edge modes is 
dependent on boundary conditions. 

To analyze the presence of edge states in gapped BLG 
we solve the wave equation for different edge structures 
(that do not couple states in different valleys), to which 
corresponding appropriate boundary conditions are asso- 
ciated. Possible boundary conditions for the single- valley 
models can be derived from the constraint of vanishing 
probability current across the boundary (we comment on 
their physical realizations shortly later). Again we use 
one single valley Hamiltonian, described by Eq. ([2]), as 
our example. Let us write the two-component wavefunc- 
tion as ^'(r) — {ipi{r),tl>2{r))'^ . From the single valley 
continuity equation ^(^^5") + V • j = 0, we obtain the 
probability current density 



A1 o B1 • A2 o B2 • 



j{r) = [^{^p'^^-'^p2 + V'2^+'0i) + c.c.]x 
+ [iipld_'ip2 - tpld+i^i) + c.c.]y , 



(15) 



where d± = dx^idy and c.c. stands for complex conjuga- 
tion. Without losing generality, let us suppose we have a 
semi-infinite sample occupying the upper half plane with 
its edge along y — Q. The requirement that no current 
should fiow through the boundary implies 

3y{v = 0) = - r2d+^Ji) + c.c.]\y^o = . (16) 

This equation can be satisfied in various settings, each 
accounting for one possible boundary condition. For ex- 
ample, we can easily see that any of the following condi- 
tions satisfy the above equation. 



^1 


y=0 = {d+^l}l)\y^O 


= ; 


(17) 


-02 


y=a = (9-V-'2)|y=0 


= ; 


(18) 




V'lla=0 = V'2|y=0 


= ; 


(19) 




y=0 = (5-V'2)|y=0 


= . 


(20) 




FIG. 4. Examples of defect-free edges (left panels) of BLG 
and their corresponding electronic spectra (right panels; bulk 
states in blue, subgap states at opposite edges in red and 
green, dispersionless edge states in gray). None of the edges 
shown here couples valleys, which ensures the single-valley 
models being valid and definite. The numbers of the subgap 
edge modes nevertheless are different for different edge config- 
urations, confirming the absence of bulk-edge correspondence 
for individual valleys. Note that in the right panels, (dis- 
persion relation) lines extending to the left and to the right 
actually connect at the Brillouin zone boundary. 



In principle, the boundary conditions obtained this 
way should be verified physically by finding correspond- 
ing terminations of the BLG lattice, which is in general a 
formidable job. As a matter of fact, the examples given 
m Eqs. pTpOl can be derived alternatively from the 
tight-binding models by considering specific sublattices 
that are missing at a sample edge. This latter approach 
follows in spirit Brey and FertigS^ who dealt with single- 
layer graphene nanoribbons, and has been discussed in 
details for the case of BLG in Rei^. At this point, 
without entering details, we simply list the physical real- 
izations of the boundary conditions ( ITpO ) in the same 
order in the left panels of Fig. |4] Note in particular that 
all the edges shown are parallel to a zigzag configuration, 
which leaves the valleys uncoupled and hence allows the 
applications of the single- valley models. 

Also illustrated in Fig. |4] are the bulk band profiles 
(blue areas) and the dispersion relations of edge states 
(red and green curves for two opposite edges) correspond- 
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ing to the specific type of edges shown in the same row of 
the figure. These spectra of electronic states can be cal- 
culated either from solving the single- valley models with 
boundary conditions discussed above, or from exact di- 
agonalization of tight-binding models including edges - 
the results of the two approaches agree quantitativeljP^. 
One can immediately see that, as a most important fact, 
the number of subgap edge modes per valley per edge ex- 
plicitly depends on the edge configurations, and can be 
1, 2 or 0. Therefore, a naive extension of bulk-edge corre- 
spondence to the individual valleys, which would predict 
the above number to be 1 independent of the specific 
boundary, does not hold for gapped BLG. 

We close this section with two remarks. First, our 
focus in this paper centers on the subgap edge states, 
which are dispersive and may contribute significantly to 
the subgap conductance. Nevertheless, there can be an- 
other set of edge states (also depending on boundary 
conditions) that are dispersionless and form flat bands^° 
(shown as gray lines in Fig. |4|. Indeed, the presence of 
both types of edge states as a whole has a well-defined 
topological origin from the bulk^S^^ and can be formulated 
in terms of the Zak's geometric phase associated with the 
Bloch states for which an explicit account of the crystal- 
lographic periodicity of the BLG lattice is necessarjE^. 
Second, it is an interesting problem to understand the 
absence of valley-specific bulk-edge correspondence in the 
context of Laughlin's gedanken experiment - the inequal- 
ity (for each valley) of the quantized Hall conductivity 
and the number of gapless edge modes seems to be a para- 
dox. The key to resolve this puzzle is to notice that the 
valley quantum number, unlike charge, is not conserved 
during the adiabatic insertion of flux which is normally 
considered in a Laughlin's gedanken experiment. Such 
non- conservation of the valley quantum number prevents 
a simple relation between the quantized valley-specific 
Hall conductivi ty and the number of valley-polarized gap- 
less edge moded^iissj 



IV. MARGINAL TOPOLOGICAL CHARACTER 
OF GRAPHENE SYSTEMS 

Although the analogy between valley-based graphene 
systems and spin-based topological insulators has sig- 
nificant limitations, the bulk topological properties of 
graphene systems, as we show below, play nevertheless 
an important role, with experimental consequences. To 
emphasize the connection of graphene and topological 
insulators, and at the same time to stress that there 
are important limitations to this connection, we refer to 
graphene as to a "marginal topological insulator". We 
borrow the term "marginal" from Volovik to character- 
ize the occurrence of noncompact momentum space in 
the case of massive Dirac fermion^^H the following 
section we illustrate that such a marginal character has 
direct experimental implications. 




FIG. 5. Illustration of the mapping g(fc) for gapped BLG. 
The left panel shows the solid angle that g covers when 
k runs over the whole plane, the double arrow signifying the 
double covering of the (upper) hemisphere. The right panel 
shows how two marginal topological mappings can be "glued" 
to form a well-defined topological mapping, which is the case 
for a BLG-BLG domain wall with opposite gap (A) signs on 
the two sides (the sub-/super-scripts + and — denote the do- 
mains), ck in the left panel case is not a topological invariant, 
but 5ck in the right panel case is. 



A. Domain wall 

The noncompact momentum space for single-valley 
models can be sometimes "cured" to obtain well- 
defined nontrivial topological properties. One such 
cas e is a smo oth topological defect, such as a domain 

^al ]ll | 12 | a9 | 4U | 

In the case of gapped BLG, a domain wall can be real- 
ized controUably by gating two adjacent areas of a sam- 
ple with opposite polarities (i.e. opposite signs of A)(^. 
The resulting nontrivial topological character that can 
be ascribed to an individual valley can be understood as 
follows. We have discussed in the previous section that 
the g vectors -corresponding to the pseudo-spin polar- 
ization of Bloch states- cover a hemisphere twice when 
k runs over the 2D fc-plane in the case of a single valley 
of uniformly-gapped BLG (see the left panel of Fig. |5| . 
It is also clear that, for a specific valley, the hemisphere 
which g covers is determined by the sign of A. It fol- 
lows that for a domain wall across which A changes sign, 
the bulk g vectors of the two domains cover two oppo- 
site hemispheres (each twice) which seamlessly connect 
on the equator at A; — )> cxd (where A in the Hamiltonians 
becomes irrelevant; see the right panel of Fig. [5]). In 
this way the momentum space -within the same valley- 
is compactified as a 2D sphere and the topological in- 
variant characterizing the mapping from such a sphere 
to the parametric g-sphere is evidently nontrivial. In 
fact, this topological invariant is given by 5cr = c^^ — c:^ 
(t = K,K'), which is the difference of the number Ci-'s 
for the two domains. 

For a BLG domain wall, therefore, 6ck/k' — ±2, 
and the difference is a topological invariant. More- 
over, it is this difference that defines the number of the 
valley-polarized chiral modes via the spectral asymmetry 
theorem^^. While the theorem readily applies to smooth 
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FIG. 6. The localization length \ioc of low-energy edge states 
as a function of edge roughness depth d (divided by ao, the 
intralayer nearest-neighbor distance of graphene). The simu- 
lations to extract Aioc are carried out starting from different 
ordered edges, with or without chemical disorder. In all cases, 
sufficiently strong disorder leads to the same value of Xioc re- 
gardless of any detail. This universal localization length, A*, 
is around 20 nm for the current case shown here, correspond- 
ing to a gap of approximately 100 meV. 



- quasiclassical - interfaces such as domain walls, its con- 
ditions are in general violated at the sharp interface be- 
tween graphene and vacuum, even when the interface 
does not mix valleys. The problem is that one cannot 
adiabatically interpolate between the spectrum of gapped 
BLG to the spectrum of vacuum, while one easily can in 
the case of an interface between the regions with A > 
and A < 0. Indeed this is what has been analytically and 
numerically verified by Martin, Blanter and MorpurgcJ^i]. 
Moreover, the presence of these "interface" gapless modes 
-unlike those discussed previously for plain BLG edges- 
does not depend on microscopic details such as the crys- 
tallographic configuration of the domain wall, which con- 
firms the topological origin of these modes. 



B. Gapped BLG with Disordered Edges 

The marginality of single- valley models implies, on the 
one hand, that a small perturbation in the Hamiltonian 
or boundary conditions can have a substantial effect on 
presence or absence of the gapless modes (e.g. see Ref.f^^ 
- this is merely a restatement of the absence of bulk-edge 
correspondence for individual valleys. On the other hand, 
considering a realistic graphene device where the edges 
of the sample are inevitably under the influence of dis- 
order, physical properties will necessarily be subject to 
ensemble averages, such that, at sufficiently strong disor- 
der, universal behaviors are expected to occur. It is thus 
interesting to ask in this context whether, or how, the 
marginal bulk property can be manifested in the subgap 
regime. 




FIG. 7. Local density of states at zero energy (center of the 
gap) in a typical two terminal device of gapped BLG with 
disordered edges. Edge states are universally present in the 
subgap regime with long localization length. 



Such a question has been addressed by Li et ahl^Slwith 
the motivation to investigate the "extra" contribution by 
the edge states to the subgap conduction in realistic ex- 
periments on gapped BLG. It is found that indeed strong 
edge disorder leads to universal presence of subgap edge 
states which are localized with fairly long localization 
length. To be specific, Fig. [6] shows how the localization 
length Aioc -extracted numerically with the aid of the 
scattering theorjSUilL depends on the strength of edge 
disorder (d denotes edge roughness depth; solid lines rep- 
resent the effect of structural disorder only and broken 
lines represent the effect of structural plus chemical disor- 
der; see Ref.^for details), with various choices of initially 
ordered (e.g. zigzag, armchair) edges, and a fixed gap 
size at an experimentally accessible value. For the cases 
where the starting edges accommodate pre-existing edge 
states (counter-propagating for two valleys), gradually 
turning on disorder leads to localization and a shrinking 
in the transverse size of the edge states (see the solid lines 
in blue and green) - this is certainly expected because 
of the lack of protection for the edge states from back- 
scattering in the system. What is unexpected is that 
when starting from initial edges where no edge state ex- 
ists even in the absence of disorder, disorder starts to in- 
troduce (localized) edge states inside the gap. Moreover 
at relatively weak disorder, increasing disorder enhances 
their Ajoc (see the solid line in red). Most interestingly, 
the localization length in all cases converges to the same 
value at sufficiently strong disorder (see the solid lines in 
the large-cZ limit, and also the broken lines where dom- 
inant chemical disorder leads to very fast convergence 
even at small d). The universal value A* depends only on 
the gap size (A* oc l/-\/A; it is around 20 nm in the cases 
illustrated in Fig. [6] corresponding to a gap of approxi- 
mately 100 meV). These values of the localization length 
can have experimental consequences in real device^^H, 

Thus subgap edge states with a long localization length 
are generically present at disordered edges of gapped 
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BLG (see Fig. [7|. These states can manifest them- 
selves in transport experiments at finite temperature. 
In particular, when the subgap transport is dominated 
by one-dimensional variable range hopping along the 
edges (which is the case if the hopping through the 
bulk disorder are sufficiently suppressed), the temper- 
ature dependence of the conductance is of the form 
G cx exp[— (r*/r)^/^], with T* the characteristic energy 
scale that determines the low temperature transport be- 
havior. For practically relevant parameters, T* is pro- 
portional to, but approximately one order of magnitude 
smaller than the size of the bulk gapl^. Indeed in very 
clean gapped BLG devices, in which the bilayer is sus- 
pended and not in contact with a substrate, an excess 
conductance, thermally activated with a characteristic 
energy one order of magnitude smaller than the gap, has 
been observed experimentalljE'. These observations are 
in agreement with our predictions. 

To verify the link between presence of the weakly local- 
ized subgap edge states and the marginal bulk properties, 
we performed a comparative study between the current 
BLG system and a half-filled square lattice with different 
on-site energies on nearest neighbors. Such a square lat- 
tice is fully-gapped but has a topologically trivial band 
structure's. It turns out that in the latter system no 
subgap edge states are present neither in the ordered nor 
in the disordered limit; subgap transport is solely due 
to weak direct tunneling between contacts which is even 
further suppressed in the presence of edge disorder. The 
sharp contrast between the two systems strongly suggests 
that the universal features in terms of the (localized) sub- 
gap edge states in gapped BLG is indeed a manifestation 
of its non-trivial marginal topological bulk properties. 



V. CONCLUSION 

To conclude, we have carefully explored the apparent 
analogy between gapped graphene system and topologi- 
cal insulators, with the central element of the correspon- 
dence being the valley degree of freedom in the former 
and electron spin in the latter. We find that the strong 
formal resemblance is fundamentally limited by the dif- 
ferent topological structures in the electron spectrum un- 
derlying the two systems. We have used gapped BLG 
as a concrete example to facilitate the comparison, and 
more importantly, to demonstrate how the "marginal" 
topological properties associated with valleys can man- 
ifest themselves significantly in experiments. Yet, such 
marginal characteristics are generic for graphene systems, 
and can be readily seen by considering other graphene 
multi- and single layers. 
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